Stochastic deformation of a thermodynamic symplectic structure 



P.O. KazinskiQ 

Physics Faculty, Tomsk State University, Tomsk, 634050 Russia 
(Dated: January 9, 2009) 

Abstract 

A stochastic deformation of a thermodynamic symplectic structure is studied. The stochastic deformation 
procedure is analogous to the deformation of an algebra of observables like deformation quantization, but 
for an imaginary deformation parameter (the Planck constant). Gauge symmetries of thermodynamics and 
corresponding stochastic mechanics, which describes fluctuations of a thermodynamic system, are revealed 
and gauge fields are introduced. A physical interpretation to the gauge transformations and gauge fields is 
given. An application of the formalism to a description of systems with distributed parameters in a local 
thermodynamic equilibrium is considered. 
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I. INTRODUCTION 



A geometric description of thermodynamic systems has a long history and was initiated by 
Gibbs |jl|. In this approach, equations of state of a thermodynamic system are represented by a 
surface in a space of thermodynamic parameters. Later this geometric formahsm was developed 
in the works [§, ^, ^, ^, 0] on axiomatic foundations of thermodynamics, where its laws were 
formulated in terms of differential forms. Due to Hermann Q symplectic and contact geometries 
acquire a distinctive form in thermodynamics, though some elements of these geometries were 
introduced by Gibbs as well. Notions of symplectic, contact, Riemannian and Finslerian geometries 
in thermodynamics get a further development in the papers 12, 13, 14, 15, 17|. Since 

the symplectic structure arose, a strong analogy of thermodynamics with classical mechanics and 



optics was understood soon afterwards |18, 19, 21, 22, 23, |25[|. Thermodynamics was realized 



as a Hamiltonian dynamical system. The next logical step is to "quantize" such a dynamical system 



in order to describe its fluctuations, as it was posed in |19]. At the same time this "quantization" 
is not a quantal one, which is given in [^], but has to result in a Fokker-Planck type equation for 
a probability distribution of thermodynamic variables. The inverse procedure of "dequantization" 
(the weak noise limit) of the Fokker-Planck equation is of a common knowledge |^, |2^, ^ and 
also endows thermodynamics with the symplectic structure. In this paper, we shall see that the 
desired "quantization" is a stochastic deformation procedure proposed in [^], which is an analog 
of the algebraic approach to quantization known as deformation quantization 



31 



A theory of fluctuations of thermodynamic quantities is well-known both for equilibrium |34] 
and non-equilibrium processes, and has a huge literature (see, e.g., [27, 35, 37, p8[ 



43, 45, ^). Thus the aim of "quantization" lies in a construction of stochastic mechanics^ by 
means of a procedure similar to an ordinary quantization, while it should give rise to the standard 
theory of fluctuations and be equivalent to operator methods of stochastic mechanics pSl 4S, 



50, 51, 52|. Stochastic deformation applied to thermodynamics complies with these requirements 
and reproduces known operator methods in some particular gauges. Gauge transformations and 
gauge fields are necessary ingredients of stochastic deformation. Furthermore, as we shall see, they 
are intrinsic to thermodynamics. Transformations of this kind arise occasionally in the literature 



| 28| , |5^ as some tricks to prove, for example, an equivalence of the Doi [U, |5g, |5^ and Martin 



^ To avoid misunderstanding we point out that we use the term "stochastic mechanics" in a general sense hke a 
notion that unifies various physical models with dynamics obeying some kind of the master equation. It is not 
Nelson's Stochastic mechanics mm. 
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Siggia-Rose ]48, 58, 30, 62, 63, 64 1 formalisms |65], or to establish certain properties of 
operator's spectrum of the master equation, or they appear in the form of nonstandard inner 



products or as the a degree of freedom in Umezawa's thermo field dynamics |67] 

and so on. As long as thermodynamics and stochastic mechanics have such symmetries, it is of 
importance to represent their governing equations in an explicitly invariant form like, for example, 
to represent the Maxwell equations in an explicitly Lorentz and gauge covariant fashion. 

Stochastic deformation of the thermodynamic symplectic structure reveals, firstly, a perfect 
analogy of thermodynamics and stochastic mechanics on the one hand with classical mechanics 
and its quantum deformation on the other hand. Secondly, it allows us to discover natural gauge 
structures of thermodynamics and stochastic mechanics and formulate them in gauge invariant 
forms. In comparison with the standard operator approach to stochastic mechanics, we effectively 
introduce an additional Stiickelberg field and corresponding gauge fields, which result in gauge in- 
variant dynamics. A significant feature of these rather formal manipulations is that the Stiickelberg 
field possesses a physical interpretation. It is an entropy of the thermodynamic system, while the 
gauge fields are external thermodynamic forces acting on it. 

It is worthy to mention the works that are closely related to the subject matter of the paper. 
A probabilistic stochastic deformation of a linear symplectic structure was studied in the papers 



within the framework of Zambrini's Euclidean quantum mechanics [39, [70[|. The latter stems 



from Schrodinger's works [71| made in attempt to give a stochastic interpretation to quantum 
mechanics. In the papers [^2|, Lagrangian and Hamiltonian descriptions of Smoluchowski diffusion 
processes are given as a special case of the general formalism based on diffusion-type equations. 
Notice also that something similar to the gauge fields we are about to consider arose in |41]. Gauge 



transformations, but of the different type, was introduced into thermodynamics in |14]. 

The paper is organized as follows. We start with a formulation of thermodynamics in terms 
of symplectic geometry (Sec. ||). Then we realize non-equilibrium thermodynamics as a certain 
Hamiltonian dynamical system and establish its gauge invariance. Here we also introduce gauge 



fields and provide their physical interpretation. In Sec. |I| we consider a stochastic deformation 
of the obtained Hamiltonian dynamical system and investigate some its properties. We introduce 



an operator of dissipation, which is the stochastic analog of a dissipation function 40, 73], and 
express an entropy production through it. As an example we apply the formalism to a thermo- 
dynamic system being in a local thermodynamic equilibrium. In conclusion, we outline prospects 
for further generalizations and research. We assume Einstein's summation rule unless otherwise 
stated. Latin indices run from 1 to d, where d is a number of extensive variables, and zero index 
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corresponds to time. 



II. THERMODYNAMICS 

All the statistical properties of a quantum system in a thermostat are determined by the par- 
tition function 

Z{pi,...,Pd-i,x'^) := Spexp ^-^Pi^'^ , (1) 

where are quantum operators of additive integrals of motion, x'^ is some fixed extensive variable, 
for instance, the volume, and pi are intensive parameters thermodynamically conjugate to x'^ or 
thermodynamic forces. The Hamiltonian and the reciprocal temperature are among these conjugate 
pairs. Differentiating the partition function with respect to pi and making the Legendre transform 
we arrive at the first law of thermodynamics 

dS{x)=pidx\ (2) 

where are averages of the corresponding operators, S{x) is the entropy of the system, which 
is the Legendre transform of the Massieu function [^, ^4|, |7^ $ := —InZ, and we introduce the 
intensive parameter p^ := ddS{x) conjugate to x'^. In geometric terms, the first law (Q) endows 
the space of states {x,p) of a thermodynamic system with a symplectic structure specified by the 
symplectic potential 9 := pjdx* — dS. According to @ the system is confined to the Lagrangian 
surface of the symplectic 2-form dO. Inversely, any Lagrangian surface uniquely projectable to the 
space of extensive variables {x*} can be locally represented by the equality of the form (||) (see, e.g., 
iPUl ). For reasonable physical systems the Lagrangian surface is uniquely projectible, otherwise 
a change of thermodynamic forces does not vary extensive parameters, that is, we have a system 
with zero (generalized) compressibility. 

In nonequilibrium with the thermostat, the thermodynamic system moves along the Lagrangian 
surface. Besides, if we directly (not by means of intensive parameters) change the entropy of 
the system, the Lagrangian surface (||) also evolves. A natural generalization of the first law to 
nonequilibrium processes looks like 

dS{t, x) = Pidx^ — H{x,p, t)dt, (3) 

where H{x,p,t) is a thermodynamic force conjugate to time or the Hamilton function^. The 

^ Of course, this Hamilton function is not related, at least immediately, to the Hamiltonian entering the partition 
function (|l|). In ksl this thermodynamic force is called the kinetic potential. 
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Hamilton function taken on the Lagrangian surface is a source of the entropy and it vanishes if the 
functional form of the entropy does not change. The non-stationary first law (^) represents the 
Hamilton-Jacobi equation. It is valid when it is reasonable to attribute certain values of intensive 
parameters to the whole thermodynamic system. 

The requirement that the system should be confined to the Lagrangian surface restricts the form 
of the Hamilton function. Here we consider two equivalent mechanisms how to keep the system on 
the Lagrangian surface 

Ti := Pi - d^S = 0. (4) 

The first one is a "rigid" method. We demand that Tj are integrals of motion of the Hamilton 
system associated with (^). A general form of a regular in momenta Hamilton function meeting 
this requirement is 

H = c{t) - dtS + Tiv'{t, x) + h'ig'^it, x)Tj + ..., (5) 

where and g^^ = g^'^ are some contravariant tensors, dots denote the terms of a higher order in 
Tj, and c{t) is some function. Redefining the entropy we eliminate the latter function. Then for 
a stationary entropy the Hamilton function (^) is zero on the Lagrangian surface. The Hamilton 
equations of motion of the thermodynamic system become 

x' = {x\ H} ^ v\ ti = dtT, + {Ti, H} ^ 0, (6) 

where curly brackets denote the Poisson brackets, and approximate equalities mean that we take 
equations on the Lagrangian surface (Q). Now we see that the vector field describes a drift of 
the system. The tensor g^^ and higher terms in the expansion (^) gain a physical meaning only 
upon stochastic deformation of the Hamiltonian system, which will be considered below. They are 
responsible for a probability distribution law of fluctuations of statistical averages. 

The Hamiltonian action functional associated with @ is obvious. If the expansion of the 
Hamilton function (P) terminates at the term quadratic in Tj and the tensor g'^^ is nondegenerate, 
it is not difficult to obtain a Lagrangian form of this action 



S[xit)] = / dt 



1 dS 

-9^j{x' -V'){X^ -V^) + — 



(7) 



where gij is the inverse of g'^K This functional is known as the Onsager-Machlup action |38, 75, 76]. 
It measures an entropy change along a trajectory x{t). The first term measures an additional 
entropy production caused by fluctuations deviating intensive parameters of the system from the 



Lagrangian surface. This term vanishes in the thermodynamic hmit @. The second term is a 
change of the thermodynamic entropy. The principle of least action corresponding to (Q) says that 
the system moves to equilibrium with a minimum of the entropy production. Therefore the first 
term has to be nonnegative for an arbitrary trajectory and the tensor g^^ must be positive definite. 



Further, if the fluctuations are described by a general Markov process, it can be proven |43] that 
dp H should be positive definite. 

Under the above restrictions on g'^^ any drift vector field can be represented in the form 

v' = g'\djS-Aj), (8) 

where Ai{t,x) is some local 1-form, which we shall call the gauge field. The quantity in the 
brackets is the total thermodynamic force, while the gauge field Ai is the external force exerting 
on the system. In particular, the Onsager principle [^] postulates the drift dH) and that the 
thermostat acts on the nonequilibrium system as 

Ai{t,x) =diS{t,xo), (9) 

where Xq is the solution of (Q) at fixed intensive parameters pi characterizing the thermostat. To 
provide a stability of the state xq of the system we have to require a negative definiteness of the 
Hessian dijS{xo), otherwise the system is in a phase transition state. If the total thermodynamic 
force is given, the relation (|8|) can be only a linear approximation to a nonlinear response of the 
system to the uncompensated force. Besides, there are systems that do not obey the Onsager 



principle (see, for physical examples, |27, 36, |6l|, ^ [78|| ) 



The first law (^) with the Hamilton function (^) is invariant under the following gauge trans- 
formations 

Pi^Pi + di^{t,x), Ai^ Ai + diC{t,x), Sit,x) ^ S{t,x)+^{t,x). (10) 

Their existence reflects the fact that a gradient part of the external force can be attributed to the 
system itself redefining its entropy. In other words, these transformations relate equivalent ther- 
modynamic systems (system)+(thermostat). They are not distinguishable within thermodynamics 
since the total forces do not change under the gauge transformations ([To|). 

Now we are in position to introduce the second method to confine the system to the Lagrangian 
surface. This method is, of course, equivalent to the first one and based on using an auxiliary 
compensating field. We postulate that the first law (^) is invariant under the gauge transformations 
([To|). In addition, we seek the Hamilton function, which does not depend on the entropy. Then a 



quadratic in momenta Hamilton function^ providing gauge invariance to (^) can be cast into the 
form of the Hamilton function of a charged particle in an electromagnetic field 



H = -Vig'^Vj - Ao =: K{p, x) - Aq, Vi := Pi - A, (11) 



where Ao(t,x) is the auxiliary field transforming under the gauge transformations as 

Ao^ Ao + dtat,x). (12) 

It keeps an invariance of ^ with respect to nonstationary gauge transformations. To confine the 
system to the Lagrangian surface determined by a given entropy function S{t, x), we have to choose 
the compensating field Aq so that the Hamilton-Jacobi equation is fulfilled: 

^0 - dtS = ^g'^idiS - A,)idjS - A,). (13) 

If we substitute from the Hamilton-Jacobi equation to the Hamilton function (pA]), we re- 
vert to (^) with the drift vector field (^) establishing the equivalence of two approaches. The 
gauge invariant total force Aq — dtS is the kinetic part K of the Hamilton function taken on the 
Lagrangian surface. These two functions, Aq — dfS and K, are different representations of the 
dissipative function ^'(X, X) introduced in |^^. The entropy production S in the whole system 
(thermostat) -|- (system) is standardly expressed in terms of the dissipation function 

dK 

S ^ V.J- = 2K, (14) 
dpi 

where we assume that the rate of an entropy change of the thermostat is —v'^Ai. 

Given K and A^ completely define the system and, in particular, its thermodynamic entropy 



through the Hamilton-Jacobi equation. From mechanical viewpoint a thermodynamic system (11) 
tends to an unstable equilibrium of the potential —Aq. That is, for any given x we take such 
an initial momentum p that the "particle" hits precisely the equilibrium point with zero velocity. 
Although this point is unstable equilibrium, it is of course an attractor on the Lagrangian surface. 
Not any mechanical Hamilton function ( |Tl[) having a nonstable equilibrium point in its potential 
can serve as the Hamilton function for some thermodynamic system even in a neighbourhood 
of this point. The "magnetic" field can freeze the particle (like in a magnetized plasma) so it 
never reaches the equilibrium. Linearizing particle's equations of motion in a small vicinity of the 



^ A generalization of these considerations to Hamilton functions of an arbitrary order in momenta is straightforward, 
but has no such a suggestive mechanical analogy. Quadratic in momenta Hamilton functions correspond to a 
Gaussian distribution law of fluctuations (5-correlated in time. 
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stationary point xq and assuming that the Hamilton function is stationary and the fields entering 
it are smooth enough, one finds that for the isotropic potential dijAQ{xo) = S,9ij{xo) the system 
never hits the point xq starting in an arbitrary point of its small vicinity only if 

\Xi\ > 4e, (15) 

for some i, where Aj are the characteristic numbers 

det(Fifc/'Fzj - Xg,j) = 0, Fij := ^^^j (xq). (16) 

In other words, the condition (15) reminds that the gauge fields Ai should enter the potential 



according to the Hamilton- Jacobi equation (|1^). 

III. FLUCTUATIONS 

Now we turn to fluctuations of statistical averages. As it was shown in [^], these fluctuations 
can be obtained by stochastic deformation of the corresponding Poisson structure. In our case, we 
shall deform a canonical symplectic structure associated with the thermodynamic system. 

Let us briefly recall some basic features of an algebraic stochastic deformation. For more details, 



an interested reader can consult the paper |2£] and the classical works on deformation quantization 
1 30, 31, 32, 33]. A commutative associative algebra of classical observables is constituted by real 
smooth functions over the symplectic space. We deform this algebra in a manner of deformation 
quantization, but with an imaginary deformation parameter, such that 

[x\pj]=u6i, (17) 

where is the real positive deformation parameter. Hats signify elements of the deformed associa- 
tive algebra and we imply the Weyl-Moyal star product fTT 



^^ = Y.-^ • • • ^""'"5a....a„F(z)a,,...,„G(z), (18) 



n! V2 

n=0 



where z = {x,p), an,bn = l,2d, the functions F{z) and G{z) are the Weyl symbols of the corre- 
sponding elements of the deformed algebra, Lo"'' is the inverse to the symplectic 2- form ojab- The 
generators and pj of the Heisenberg-Weyl algebra correspond to extensive and intensive param- 
eters of the thermodynamic system. The deformation parameter is not the Planck constant. It 
characterizes a variance of fluctuations and, as we shall see, is equal to doubled the Boltzmann 
constant, u = 2kB, for thermal fluctuations. 
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Another necessary ingredient of the deformation procedure is the trace functional Sp, which is a 
hnear functional on the deformed algebra mapping to real numbers and vanishing on commutators. 
An explicit formula for the trace of an element F has the form 

A state of a stochastic system is characterized by an element p with a unit trace 

Spp = l. (20) 

The pure state is specified by an additional idempotency requirement 

p2 = p. (21) 

An average of some observable F over the state p is defined by the standard formula 

(F) := Sp(/5F). (22) 

Thus, for a correct probabilistic interpretation the state p should satisfy 

{5'^{x' - x')) > 0, Vx G R'^. (23) 

The dynamics of the stochastic system in the state p are generated by the element H of the de- 
formed algebra, which corresponds to the Hamilton function H{t, x,p), and obey the von Neumann 
equation 

i^f>=[H,p]. (24) 

In the case of a linear symplectic space it is useful to realize the Heisenberg-Weyl algebra as 
operators acting in the linear space V of smooth real functions on the configuration space. Then, 
in Dirac's notations, the pure state is represented by^ 

p=\i^){f\, (v^|V) = i, W^V, {^\eV\ (25) 

where the standard inner product is understood. So, the pure state is specified by two real func- 
tions on the configuration space. For such states the von Neumann equation is equivalent to two 
Schrodinger-Zambrini (SZ) equations |6^, |7l[ 

udt\iP) = H\iP), udt{^\ = -{ip\H. (26) 



Similar projectors also arise when describing projected processes 



Hereinafter, for simplicity, we restrict ourself to the case of at most quadratic in momenta Hamil- 
tonians. Besides, we take the metric tensor gij to be a constant matrix. 

After introducing the stochastic phase S{t, x) and the probabihty density function p{t, x) to 
find a system with the values of extensive parameters 

ijj = pexY){—u~^ S), (f = exp{i'~^ S), (27) 

the operators of the total forces = — have the averages 

{V^) = jd'^x^{t,x)[-i^df,-A^{t,x)]^it,x) = {d^S-A^), ^ = M, (28) 

where, for /i = 0, we have used the equations of motion (p6|). Matching thermodynamics with 
its deformation, we should identify the phase S with the thermodynamic entropy 5. Then the 
SZ equations (|2^) are invariant with respect to the gauge transformations ( p!o[ ) both with the 
Hamiltonian (^) and (^ij). 

Consider stochastic deformation of the Hamiltonian dynamical system (|5|). With the above 
mentioned identification, one of the SZ equations (|2^) is identically satisfied 

udtiifl = -{^\{-dtS + fiv\t,x) + ^f.g'^fj), (29) 

and the other equation becomes the Fokker-Planck equation describing fluctuations of the thermo- 
dynamic system 

dtp = -d,{-'^g'^djp + v'p). (30) 

We see that the coefficients of expansion (^) in terms of T, are nothing but the cumulants of the 
probability distribution of fluctuations. 

Deforming the Hamiltonian dynamics generated by (|Tl|) we arrive at two equations: the Fokker- 
Planck equation (|30|), and the Hamilton- Jacobi equation with stochastic correction or the 
Burgers equation |]70| , 

^0 - dtS = ^g'\^^S - A,){d,S - A,) + '^^^[g'^{^,S - A,)], (31) 

defining the dissipation function Aq — dtS. The form of gauge transformations ( |To| ) and the rep- 
resentation (^^ show that the gauge group is the Abelian one dimensional Lie group isomorphic 
to 5'0(1,1). The SZ equations can be rewritten in an explicitly covariant form with respect to 
the gauge transformations if we group the two functions i/j{x) and ip{x) into one vector ^'"(a;) and 
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define the pseudo-Euclidean scalar product of such vectors as 
Then the SZ equations (^) look like the matrix Schrodinger equation 



1 




tp{x) 


1 







(32) 



H 




(33) 



where the cross denotes a conjugation with respect to the standard inner product. The evolution 
is generated by the matrix Hamiltonian, which is skew-adjoint with respect to the scalar product 
^). On introducing the self-adjoint covariant derivatives 



where 



UJr- 



^l|3 ■- 



1 
-1 



T 
UJ T] 



(34) 



(35) 



the SZ equations with the Hamiltonian (^) read as 



(36) 



In this form, the SZ equations are very similar to the quantum mechanical Schrodinger equation. 
The difference is that the generator lv of the Lie algebra so(l, 1) appears instead of the generator 
of u{\) ~ so(2). Further, the equations of motion (|3^ ) are immediately generalized to a non- 
Abelian case. If we have identical thermodynamic systems in a thermostat then the global 
symmetry group for a whole system will be SO{N,N) D Sp{N). The symplectic group Sp{N) 
arises since in the case of identical systems the matrix rjuj is the unit symplectic matrix, which 
must be preserved by the symmetry transformations. A detailed investigation of peculiarities of 
non-Abelian systems and their physical interpretation will be given elsewhere. Here we just notice 
that the group SU (1,1) locally isomorphic to SO{2, 2) n Sp{2) /50(1, 1) was studied in the context 

, p^]. The pair {r],uj) is the analog of an almost generalized 



of generalized coherent states pO, 81 



product structure on the Whitney sum TM © T*M (see, for example, 

Let us consider how some standard thermodynamic relations look in our framework. The 
condition of a detailed balance in some state |?/')('/^| looks like 



ip{x) o H o il){x) = i/^ix) o o ip{x), 
11 



(37) 



where o means a composition of operators. That is, we can make the gauge transformation so 



that (/9 = -0, the Hamiltonian H = having a nonpositive spectrum (see, e.g., |2^, |3^, 51, |85|). 
Then the stochastic phase is identified with a half of the entropy of the whole system (thermo- 
stat) + (system). For Hamiltonians of the form ( |Tl| ) with nondegnerate metric (7'-' the probability 
density function p(x) in this state^ is proportional to exp[2i^^-'^S(x)], where S(x) := S{x) —piX^, 
provided the Onsager principle (^) is fulfilled. Fixing the extensive variable x'^ and applying the 
WKB-method, we obtain that the leading in u terms of the characteristic function In z{p,x'^) of 
the probability distribution p at fixed take the form 

In z{p, x'^) = 2v"^ In Z(p, x'^) + Sp In |a^^'$(p, x'^) |/2 + . . . (38) 

Thus we arrive at the well-known result that in the leading order the correlators of statistical 
averages are proportional to the correlators computed with the help of the partition function (||). 
In the Heisenberg representation, the operator of the system entropy change is 

S = dtS + u-\S,H]=dtS + diS—, (39) 

OPi 

where diSd/dpi is a differentiation of the Heisenberg- Weyl algebra acting on its generators in 
an obvious manner. Consequently, the gauge invariant entropy production in the whole system 
(thermostat) + (system) is given by 

t:={d,S-A,)^, (40) 
dpi 

where the dot is just a notation in the case of Fij 0. The relation (^) is the stochastic 
(noncommutative) analog of its thermodynamic counterpart (14). It is reasonable to define the 
operator of a purely fluctuational entropy production as 

-£'f, = -u-^[x\H]fi, (41) 

that follows from a path-integral representation of its average. This kind of the entropy production 
disappears in the thermodynamic limit. If the system with the Hamiltonian ( pT] ) is in the state, 
where the detailed balance takes place, the average entropy production S is zero. The average of 
the fluctuational term becomes 

- u~^{[x\ k]f,) = -u{g'^dijT.) = 2(5*^a,Sa,S). (42) 



* The use of this probability density function in the study of fluctuations of the Van der Waals gases can be found, 
for example, in ^ . 
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It is positive and of the order of u. This is vahd in the leading order in u for any Hamiltonian 
that allows a detailed balance. Thus the entropy production S is caused by equalizing large 
(macroscopic) differences of intensive parameters of the system and the thermostat, whereas the 
fluctuational term is responsible for the entropy change made by equalizing small (microscopic) 
differences of intensive parameters originating from fluctuations. The microscopic deviations can 
be estimated from the well-known thermodynamic uncertainty relation (see, e.g., ||T^). In the state 
with a detailed balance, we have (for a proof see, e.g., p9| ) 

{{x'f){{di^f)>v'^/A, (no summation). (43) 

A natural generalization of the above considerations to more realistic nonequilibrium systems 
possessing spatial gradients of intensive parameters is straightforward. We break the system to 
subsystems of a fixed volume, which are small enough to have homogeneous intensive parame- 
ters and sufficiently large to apply a thermodynamic description to them, i.e., the system is in 
a local thermodynamic equilibrium. Then the simplest model following from first principles of 



thermodynamics (^) and (^ prescribes a diffusion- like evolution ||^, 14, 35, 44, B8| 



dtnt,x) = j dyg^\x,y) (^^^ - ^|^) , g'-'ix^y) = 9^\y,x) = g''^{x,y), (44) 

where (j)"'{t,x) are densities of the extensive variables except the volume, the functional S[(j)\ is the 
thermodynamic entropy of the whole system, and g""^ is some positive definite matrix for any x 
and y, which measures a linear response of the extensive variable a of the subsystem in the point 
x on a difference of thermodynamic forces b of the subsystems located at x and y. The use of the 
linear response relation is justified by small, by construction, deviations of intensive parameters 
of neighboring subsystems. The total values of extensive variables are conserved by the evolution 



(44). The Hamilton functional of the form (|ll| ) corresponding to Gaussian fluctuations becomes 



/7[t,0,7r] = i j dxdyiTTaix) - TTaiy))g''\x,y){7rb{x) - TThiy)) - Ao[t,(l)], (45) 
where tTq are the intensive parameters canonically conjugate to 0". If external thermodynamic 



force fields are applied, the momenta should be replaced by the covariant derivatives (11). The 
Hamiltonian formalism is also preferred to the Lagrangian one as the Onsager-Machlup action 
is nonlocal for local g""^. The functional of a thermodynamic entropy increases with the evolution 
(44) and acquires maximum when all the intensive parameters become homogeneous. 

Upon stochastic deformation of the model ( |45| ) we see that the fluctuating system possesses the 
states in which the detailed balance (^) takes place: 

^/ifocexp J Xa(p"'{x)dx, (46) 
13 



where are some constants and the stochastic phase vhup is the thermodynamic entropy. We 
divide the whole system into a thermostat and the system in it imposing the constraints 

{^\{T:a{x)-pa) = =0. (47) 

Here x runs points of the thermostat and pa are fixed values of its intensive parameters. These 
constraints are preserved by the evolution and just say that the intensive variables of the thermostat 
do not fluctuate. The Hamiltonian corresponding to (|45| ) is Hermitian and has a nonpositive 
spectrum. Hence any state of the system satisfying ( |47| ) tends to the ground state described by 
the probability density functional of the expected form 



/?[(/)] cx exp 



2u-^ (Ssys[(t>\ - [ dxparix)) 

V J sys / 



(48) 



'sys 

where the integral is taken over the system in the thermostat and Ssys[<P] is the thermodynamic 
entropy of this system. As before, the average entropy production S relative to this state vanishes. 
For local g"''' the density of fluctuational entropy production is zero in the thermostat, while it is 
of the form (E^) in the system. As a matter of fact, the matrix g"'^{x,y) depends on the fields (/)"'. 



A generalization to this case is easily realized along lines |29, 41, w 



IV. CONCLUDING REMARKS 

Let us mention some possible modifications and generalizations of the formalism evolved here. 
Notwithstanding we distinguish extensive and intensive parameters, we did not actually use these 
properties. If the thermodynamic system possesses the "gauge" symmetry 

x'diS = S, (49) 

i.e. S{x) is a homogeneous function, a division into intensive / and extensive E variables is achieved 
by 

51 := {x'T„ 1} « 0, SE := {x%, E} ^ E. (50) 

As it follows from ( ^9|) the Hessian is degenerate. Therefore we have to fix the extensive parameter 
x'^ in (^) and work in the sector of remaining extensive variables and their conjugate. The fixed 
extensive parameter and its conjugate are expressed in terms of the independent ones by means 



of (49). To put it in another way, we fix the "gauge". The division into extensive and intensive 



variables restricts admissible canonical transformations of the phase space of a thermodynamic 
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system. Such a separation is preserved by arbitrary changes of intensive parameters and only 
hnear transformations (with coefficients depending on the intensive variables) of extensive ones. 
In geometric terms, we realize the phase space as a Lagrangian fiber bundle |2^] with the base 
parameterized by intensive variables and the fiber represented by the configuration space. 

Noteworthily, the thermodynamic (mean field) approach to phase transitions is naturally treated 
in the framework of symplectic geometry. A set of singular points of the Lagrangian surface, 
where the Hessian dijS is degenerate, projected to the base of the Lagrangian bundle is called 
in optics and mechanics the caustic. In our case, the degeneracy condition bounds a region of 
extensive parameters, where the system is unstable, metastable states being regarded as stable. 
The caustics surround phase equilibrium curves (surfaces) on the base. When the system moves 
to the phase equilibrium curve and intersects the caustic, an additional local maximum in the 
entropy S(x) = S{x) — piX^ of the whole system appears. This maximum corresponds to a new 
phase and becomes equal to the entropy maximum of the old phase on the equilibrium curve. It is 
interesting to apply developed methods of symplectic geometry to a study of topology and general 



forms of singularities M, 2C, pq, B3, BH, Bll of Lagrangian surfaces, and, consequently, caustics and 



phase transition surfaces. The singularities of Lagrangian surfaces, which are stable with respect 



to small deformations respecting the Lagrangian bundle structure, are classified in |21] up to ten 
dimensional phase spaces. The normal forms of Lagrangian surface singularities in phase spaces of 
higher dimensions have moduli. 

Sometimes it is useful to define an entropy of a thermodynamic system not as the Legendre 



transform of but as follows (see, e.g., |43|, |74| |92| , 



S{x) = In 



J dpexp ^^^PjX* — <I> 



(51) 



where contours of integration in complex planes are taken so that the integral converges. In a 
macroscopic limit, this definition coincides with the standard one, what is easy to see using the 
WKB-method. The extensive variables in ( |5l|) are not averaged over the ensemble as in (|^. 
Rather, they describe one system in it. The intensive variables characterizing this system are 
taken, by definition, to be pi := diS. They slightly differ from the intensive parameters pi of the 
thermostat: 

diS = diS + ^a^.fcSaf <!> + . . . (52) 

Assuming the Onsager principle (^ is fulfilled, the system with the Hamiltonian ( pAj ) and the en- 
tropy ( ^ll ) decays to the state, where a detailed balance takes place. In this state, the characteristic 
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function ( pSD is strictly proportional to \nZ up to an irrelevant additive constant. 

An entropy of a thermostat is usually described by the term —piX^ in the total entropy 
since the intensive parameters of the thermostat are assumed to be nonfluctuating. If such fluc- 
tuations become relevant, they can be naturally taken into account in the formalism of stochastic 
deformation by introducing mixed states. These states are sums of pure states with some weights 

p 

The function cr{p) is proportional to the entropy of the thermostat, while —piX^ mentioned above 
is the entropy of interaction, by analogy with the interaction part of an action functional. The 
diagonal element of ( |53D is proportional to the conditional probability Pp{x) = (p{x)ipp{x). Though 
mixed states are not exhausted by those ones. 

The analogy with mechanics suggests also possible generalizations of a simple Hamiltonian 
model (11). For example, it is interesting to consider Hamiltonian dynamics on a nonlinear sym- 



plectic space describing a thermodynamic system. Passing into the Darboux coordinates we see 
that a noncanonical symplectic structure results effectively in a changing of thermodynamic forces 
and the probability distribution law of fluctuations. Noncanonical symplectic structures appear 
naturally 44, 95, ^ in the case, when small deviations from the Gibbs distribution (jl|) 



exist 

x' = d;^ + eip), (54) 

where is a small nongradient vector field. This equality defines a Lagrangian surface of a 
noncanonical symplectic structure with magnetic fields. For example |^6|, if the density matrix of 
the system is proportional to exp[— pjX* — \{piX^Y /2 + 0(A^)], where A is a small parameter, then 
in the leading orders 

x' = di^^-\[d'^^Qp,-di^odl^oP3 + \df^oP3Pk-d'^^Qd^p^^^^ (55) 

where $o •= ~ lii.^o with Zq taken from (||). The last term in the brackets is nongradient. Defining 
a thermodynamic entropy as the Legendre transform of the Massieu function we arrive at the 
first law of thermodynamics with magnetic fields. The gradient part can be absorbed into the 
entropy by a gauge transformation. Another evident generalization of the model ([Tl| ) consists in 



introducing non-Abelian gauge fields. We only touched the problem in Sec. Ill and formulated the 



"matter" dynamics. The next step is to introduce the action functional for the gauge fields. The 
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gauge symmetries of gauge fields' action seem to be spontaneously broken, though a form of the 
action and a mechanism of the symmetry breaking are the subjects for further research. 

Besides, we saw that quantum and stochastic mechanics differ from each other by a symmetry 
group only. Namely, the Heisenberg-Weyl algebra is a Lie algebra with commutation relations 

[ei, 62] = 63, [ei, 63] = [e2, 63] = 0, (56) 

where ei, 62 and 63 are its generators and, for brevity, we consider a two-dimensional phase- 
space. This algebra includes an Abelian ideal spanned on 63. A general form of the Lie group 
corresponding to this ideal is [/"(I) x SO''{l, 1). In the two simplest cases n = 1, k = and 77, = 0, 
A; = 1 we obtain quantum and stochastic mechanics respectively. So, one can speculate about their 
unification by introducing a larger group containing the subgroups ^7(1) and 50(1, 1). 

Acknowledgments 

I am grateful to A.V. Shapovalov and A. A. Sharapov for reading a draft of this paper, discussions 
and useful suggestions. This work was supported by the RFBR Grant No. 09-02-00723 and Support 
of Russian Scientific Schools Grant No. SS-871.2008.2. 



[1] J.W. Gibbs, The Collected Works Vol. I: Thermodynamics (Longmans, Green and Co., New York, 
1928). 

[2] C. Caratheodory, Untersuchungen ueber die Grundlagen der Thermodynamik, Math. Ann. 67, 355 
(1909). 

[3] A. Lande, Axiomatische Begriindung der Thermodynamik durch Caratheodory (Springer- Verlag, 
Berlin, 1926). 

[4] T. Ehrenfest-Afanassjewa, Die Grundlagen der Thermodynamik (E. J. Brill, Leiden, 1956). 

[5] L. Tisza, Generalized Thermodynamics (MIT Press, Cambridge, Massachusetts, 1966). 

[6] J.M. Jauch, On a new foundation of equilibrium thermodynamics, Found. Phys. 2, 327 (1972); Ana- 

lytical thermodynamics. Part I. Thermostatics - general theory, Found. Phys. 5, 111 (1975). 
[7] H.B. Callen, Thermodynamics and Introduction to Thermostatistics (John Wiley & Sons, Singapore, 

1985). 

[8] R. Hermann, Geometry, Physics and Systems (Marcel Dekker, New York, 1973). 
[9] F. Weinhold, Metric geometry of equilibrium thermodynamics, J. Chem. Phys. 63, 2479 (1975). 
[10] P. Salamon, E. Ihrig, and R.S. Berry, A group of coordinate transformations which preserve the metric 
of Weinhold, J. Math. Phys. 24, 2515 (1983). 



17 



[11] R. Mrugala, J.D. Nulton, J.Ch. Schon, and P. Salamon, Statistical approach to the geometric structure 
of thermodynamics, Phys. Rev. A 41, 3156 (1990); Contact structure in thermodynamic theory, Rep. 
Math. Phys. 29, 109 (1991). 

[12] G. Ruppeiner, Riemannian geometry in thermodynamic fluctuation theory. Rev. Mod. Phys. 67, 605 
(1995). 

[13] L. Benayoun, Methodes geometriques pour I'etude des systemes thermodynamiques et la generation 

d'equations d'etat, Ph.D. thesis, Grenoble Institute of Technology, 1999. 
[14] R. Balian, P. Valentin, Hamiltonian structure of thermodynamics with gauge, Eur. Phys. J. B 21, 269 

(2001), cond-niat/0007292. 

[15] S. Preston, J. Vargo, Indefinite metric of R. Mrugala and the geometry of thermodynamical phase 
space, math/0509267vl. 

[16] J.L. Alvarez, H. Quevedo, and A. Sanchez, Unified geometric description of black hole thermodynamics, 
Phys. Rev. D 77, 084004 (2008). 

[17] S.I. Vacaru, Locally anisotropic kinetic processes and thermodynamics in curved spaces, Ann. Phys. 
(NY) 290, 83 (2001), arXiv:gr-qc/0001060. 

[18] M.A. Peterson, Analogy between thermodynamics and mechanics. Am. J. Phys. 47, 488 (1979). 

[19] R. Gilmore, Catastrophe Theory for Scientists and Engineers (John Wiley & Sons, New York, 1981). 

[20] T. Poston, I.N. Stewart, Catastrophe Theory and Its Applications (Pitman, London, 1978). 

[21] V.I. Arnold, Mathematical Methods of Classical Mechanics (Springer- Verlag, New York, 1989), 
Appedix; V.I. Arnold, S.M. Gusein-Zade, and A.N. Varchenko, Singularities of Differentiable Maps 
Vol. I: The Classification of Critical Points, Caustics, and Wave Fronts (Birkhauser, Boston, 1985). 

[22] M. Grmela, H.Ch. Ottinger, Dynamics and thermodynamics of complex fluids. I. Development of a 
general formalism, Phys. Rev. E 56, 6620 (1997); H.Ch. Ottinger, M. Grmela, Dynamics and ther- 
modynamics of complex fluids. II. Illustrations of a general formalism, Phys. Rev. E 56, 6633 (1997); 
M. Grmela, Reciprocity relations in thermodynamics, Physica A 309, 304 (2002); H.Ch. Ottinger, 
Nonequilibrium thermodynamics for open systems, Phys. Rev. E 73, 036126 (2006); H.Ch. Ottinger, 
Beyond Equilibrium Thermodynamics (John Wiley & Sons, Hoboken, 2005). 

[23] S. Sieniutycz, Hamilton- J acobi-Bellman framework for optimal control in multistage energy systems, 
Phys. Rep. 326, 165 (2000). 

[24] K. Gambar, F. Markus, Hamilton-Lagrange formalism of nonequilibrium thermodynamics, Phys. Rev. 
E 50, 1227 (1994); Onsager's regression and the field theory of parabolic transport processes, Physica A 
320, 193 (2003); F. Markus, K. Gambar, Generalized Hamilton- Jacobi equation for simple dissipative 
processes, Phys. Rev. E 70, 016123 (2004). 

[25] S.G. Rajeev, Quantization of contact manifolds and thermodynamics, Ann. Phys. (NY) 323, 768 (2008), 
arXiv:math-ph/0703061v2; A Hamilton- Jacobi formalism for thermodynamics, Ann. Phys. (NY) 323, 
2265 (2008), arXiv:0711.4319v2; A. Joseph, S.G. Rajeev, Hamilton- Jacobi formalism for string gas 
thermodynamics, arXiv:0807.3957vl; A. Coutant, S.G. Rajeev, Quantum thermodynamics of non-ideal 



18 



gases, arXiv:0807.4632vl. 

[26] R. Kubo, K. Matsuo, and K. Kitahara, Fluctuation and relaxation of macrovariables, J. Stat. Phys. 9, 
51 (1973). 

[27] N.G. van Kampen, Stochastic Processes in Physics and Chemistry (Elsevier, Amsterdam, 2007). 

[28] R. Graham, D. Roekaerts, and T. Tel, Integrability of Hamiltonians associated with Fokker-Planck 
equations, Phys. Rev. A 31, 3364 (1985); R. Graham, T. Tel, Nonequilbrium potential for coexisting 
attractors, Phys. Rev. A 33, 1322 (1986); Nonequilbrium potentials for local codimension-2 bifurcations 
of dissipative Bows, Phys. Rev. A 35, 1328 (1987); T. Tel, R. Graham, and G. Hu, Nonequilbrium 
potentials and their power-series expansions, Phys. Rev. A 40, 4065 (1989). 

[29] P.O. Kazinski, Fluctuations as stochastic deformation, Phys. Rev. E 77, 041119 (2008), 
arXiv:0711.3644. 

[30] F. A. Berezin, Quantization, Mathematics of the USSR-Izvestiya 8, 1109 (1974). 

[31] F. Bayen, M. Flato, C. Fronsdal, A. Lichnerowicz, and D. Sternheimer, Quantum mechanics as a defor- 
mation of classical mechanics, Lett. Math. Phys. 1, 521 (1975); Deformation theory and quantization, 
land II, Ann. Phys. Ill, 61 (1978). 

[32] B.V. Fedosov, A simple geometrical construction of deformation quantization, J. Differential Geom. 
40, 213 (1994). 

[33] M. Kontsevich, Defamation quantization of Poisson manifolds. I, Lett. Math. Phys. 66, 157 (2003), 
arXiv:q-alg/9709040. 

[34] A. Einstein, Theorie der Opaleszenz von homogenen Fliissigkeiten und Fliissigkeitsgemischen in der 

Nahe des kritischen Zustandes, Ann. Phys. (Leipzig) 33, 1275 (1910). 
[35] L. Onsager, Reciprocal relations in irreversible processes. I and II, Phys. Rev. 37, 405 (1931); 38, 2265 

(1931). 

[36] G. Nicolis, I. Prigogine, Self-Organization in Non-Equilibrium Systems (John Wiley & Sons, New York, 
1977). 

[37] S.R. de Groot, P. Mazur, Non-Equilibrium Thermodynamics (Dover, New York, 1984). 

[38] L. Onsager, S. Machlup, Fluctuations and irreversible processes, Phys. Rev. 91, 1505 (1953); S. 

Machlup, L. Onsager, Fluctuations and irreversible processes. II. Systems with kinetic energy, Phys. 

Rev. 91, 1512 (1953). 

[39] R. Kubo, M. Toda, and N. Hatshitsume, Statistical Physics II (Springer- Verlag, Berlin, 1985). 
[40] L.D. Landau, E.M. Lifshits, Statistical Physics. Part I (Pergamon, New York, 1980). 
[41] R. Graham, Covariant formulation of non-equilibrium statistical thermodynamics, Z. Phys. B 26, 397 
(1977). 

[42] G.N. Bochkov, Yu.E. Kuzovlev, Contribution to the general theory of thermal Euctuations in nonlinear 
systems, Zh. Eksp. Teor. Fiz. 72, 238 (1977) [Sov. Phys. JETP 45, 125 (1977)]; Fluctuation-dissipation 
relationships for nonequilibrium processes in open systems, Zh. Eksp. Teor. Fiz. 76, 1071 (1979) [Sov. 
Phys. JETP 49, 543 (1979)]. 



19 



R.L. Stratonovich, Nonlinear Nonequilibrium Thermodynamics (Springer, New York, 1992). 
D. Zubarev, V. Morozov, and G. Ropke, Statistical Mechanics of Nonequilibrium Processes Vol.1: 
Basic Concepts, Kinetic Theory (Akademie Verlag, Berlin, 1996); Vol.2: Relaxation and Hydrodynamic 
Processes (Akademie Verlag, Berlin, 1997). 

D. Jou, J. Casas- Vazquez, and G. Lebon, Extended irreversible thermodynamics revisited (1988-98), 
Rep. Prog. Phys. 62, 1035 (1999). 

V.L. Berdichevsky, Structure of equations of macrophysics, Phys. Rev. E 68, 066126 (2003). 

E. Nelson, Quantum Fluctuations (Princeton University Press, Princeton, New Jersey, 1985). 
M. Namiki, Stochastic Quantization (Springer- Verlag, Berlin, 1992). 

G. Parisi, Statistical Field Theory (Addison- Wesley, Menlo Park, 1988). 

H. Risken, The Fokker-Planck Equation (Springer, Berlin, Heidelberg, 1989). 

J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Claredon Press, Oxford, 1996). 
A.N. Vasil'ev, The Field Theoretic Renormalization Group in Critical Behavior Theory and Stochastic 
Dynamics (Chapman & Hall/CRC, Boca Raton, 2004). 

A. Suarez, J. Ross, B. Peng, K.L.C. Hunt, and P.M. Hunt, Thermodynamic and stochastic theory of 
nonequilibrium systems: A Lagrangian approach to fluctuations and relation to excess work, J. Chem. 
Phys. 102, 4563 (1995). 

M. Schonberg, Application of 2nd quantization methods to the classical statistical mechanics, Nuovo 
Cimento 9, 1139 (1953); 10, 419 (1953); A general theory of the 2nd quantization methods, Nuovo 
Cimento 10, 697 (1953). 

M. Doi, Second quantization representation for classical many-particle system, J. Phys. A 9, 1465 
(1976). 

D.C. Mattis, M.L. Glasser, The uses of quantum field theory in diffusion-limited reactions. Rev. Mod. 
Phys. 70, 979 (1998). 

P.T. Muzy, S.R. Salinas, A.E. Santana, and T. Tome, Mario Schonberg and the introduction of Fock 
space in classical statistical physics, Revista Brasileira de Ensino de Fisica 27, 447 (2005). 
P.C. Martin, E.D. Siggia, and H.A. Rose, Statistical dynamics of classical systems, Phys. Rev. A 8, 
423 (1973). 

R. Phythian, The operator formalism of classical statistical dynamics, J. Phys. A 8, 1423 (1975); 
Further application of the Martin, Siggia, Rose formalism, J. Phys. A 9, 269 (1976). 
H.K. Janssen, On a Lagrangean for classical field dynamics and renormalization group calculations of 
dynamical critical properties, Z. Phys. B 23, 377 (1976). 

C. De Dominicis, L. Peliti, Field-theory renormalization and critical dynamics above Tc: Helium, 
antiferromagnets, and liquid-gas systems, Phys. Rev. B 18, 353 (1978). 

L. Garrido, D. Lurie, and M. San Miguel, Stochastic quantization and detailed balance in Fokker-Planck 
dynamics, J. Stat. Phys. 21, 313 (1979). 

L.D. Adzhemyan, A.N. Vasil'ev, and Yu.M. Pis'mak, Renormalization-group approach in the theory 



20 



of turbulence: The dimensions of composite operators, Tear. Mat. Fiz. 57, 268 (1983) [Theor. Math. 
Phys. 57, 1131 (1983)]. 

[64] A.N. Jordan, E.V. Sukhorukov, and S. Pilgram, Fluctuation statistics in networks: A stochastic path 
integral approach, J. Math. Phys. 45, 4386 (2004). 

[65] A. Lefevre, G. BiroU, Dynamics of interacting particle systems: stochastic process and field theory, J. 
Stat. Mech.: Theory Exp. (2007) P07024. 

[66] K. Kawasaki, in Phase Transitions and Critical Phemomena Vol. 2, edited by C. Domb, M.S. Green 
(Academic, London, 1972). 

[67] T.S. Evans, I. Hardman, H. Umezawa, and Y. Yamanaka, Heisenberg and interaction representations 
in thermo field dynamics, J. Math. Phys. 33, 370 (1992); P. A. Henning, H. Umezawa, Diagonalization 
of propagators in thermo field dynamics for relativistic quantum fields, Nucl. Phys. B 417, 463, 1994; 
H. Umezawa, Advanced Field Theory: Micro, Macro and Thermal Physics (AIP, New York, 1993). 

[68] J.C. Zambrini, Feynman integrals, diffusion processes and quantum symplectic two-forms, J. Korean 
Math. Soc. 38, 385 (2001); P. Lescot, J.C. Zambrini, Probabilistic deformation of contact geometry, 
diffusion processes and their quadratures, Seminar on Stochastic Analysis, Random Fields and Appli- 
cations V, 203, Progr. Probab., 59, Birkhauser, Basel, 2008. 

[69] J.C. Zambrini, Stochastic mechanics to E. Schrddinger, Phys. Rev. A 33, 1532 (1986); Euclidean 
quantum mechanics, Phys. Rev. A 35, 3631 (1987). 

[70] Ph. Blanchard, P. Garbaczewski, Natural boundaries for the Smoluchowski equation and afElated diffu- 
sion processes, Phys. Rev. E 49, 3815 (1994); P. Garbaczewski, R. Olkiewicz, Why quantum dynamics 
can be formulated as a Markov process, Phys. Rev. A 51, 3445 (1995); P. Garbaczewski, G. Kondrat, 
and R. Olkiewicz, Burgers' flows as Markovian diffusion processes, Phys. Rev. E 55, 1401 (1997). 

[71] E. Schrodinger, Sur la theorie relativiste de I'electron et I'interpretation de la mecanique quantique, 
Ann. Inst. Henri Poincare 2, 269 (1932). 

[72] P. Garbaczewski, Differential entropy and dynamics of uncertainty, J. Stat. Phys. 123, 315 (2006); 
Modular Schrddinger equation and dynamical duality, Phys. Rev. E 78, 031101 (2008), arXiv:0805.1536. 

[73] J.W. Strutt (Lord Rayleigh), Some general theorems relating to vibrations, Proc. Math. Soc. London 
4, 357 (1873); The Theory of Sound Vol. I (Dover, New York, 1945). 

[74] R. Balian, From Microphysics to Macrophysics Vol. I: Methods and Applications of Statistical Physics 
(Springer- Verlag, Berlin, 1991). 

[75] L.M. Martyushev, V.D. Seleznev, Maximum entropy production principle in physics, chemistry and 
biology, Phys. Rep. 426, 1 (2006). 

[76] Y. Takahashi, S. Watanabe, The probability functionals (Onsager-Machlup functions) of diffusion pro- 
cesses, Lect. Notes in Math. 851, 433 (1981). 

[77] W. Ebeling, Strukturbildung Bei Irreversiblen Prozessen (Teunber, Leipzig, 1976). 

[78] K. Turitsyn, M. Chertkov, V.Y. Chernyak, and A. Puliafito, Statistics of entropy production in lin- 
earized stochastic systems, Phys. Rev. Lett. 98, 180603 (2007). 



21 



[79] F.A. Berezin, Method of Second Quantization, (Academic Press, New York, London, 1966); Feynman 

path integrals in a phase space, Sov. Phys. Usp. 23, 763 (1980). 
[80] A.M. Perelomov, Generalized Coherent States and Their Applications (Springer- Verlag, Berlin, 1986). 
[81] M. Rasetti, Coherent states and partition function, Int. J. Theor. Phys. 14, 1 (1975). 
[82] Sh. Luo, SU(1,1) coherent states and associated Wick symbol calculus, J. Math. Phys. 38, 3478 (1997). 
[83] X. Wang, B.C. Sanders, and Sh. Pan, Entangled coherent states for systems with SU(2) and SU(1,1) 

symmetries, J. Phys. A 33, 7451 (2000). 
[84] M. Zabzine, Lectures on generalized complex geometry and super symmetry, Archivum Mathematicum 

(Supplement) 42, 119 (2006), arXiv:hep-th/0605148. 
[85] H. Haken, Sinergetics (Springer- Verlag, Berlin, 1977). 

[86] Yu.Y. Beletsky, A.I. Bugrij, and A. A. Trushevsky, Phase transitions and fluctuations in ultrarelativistic 
matter with the Van der WaaJs equation of state, Z. Phys. C 10, 317 (1981). 

[87] N.V. Agudov, A.N. Malakhov, The kinetics of liquid-gas phase transitions of a Van der WaaJs substance 
with fluctuations taken into account. Chaos 4, 665 (1994). 

[88] M. Le Bellac, F. Mortessagne, and G.G. Batrouni, Equilibrium and Non-Equilibrium Statistical Ther- 
modynamics (Cambridge University Press, New York, 2004). 

[89] R. Thom, Structural Stability and Morphogenesis (Benjamin- Addison Wesley, New York, 1975). 

[90] S. Janeczko, A note on singular Lagrangian submanifolds, Ann. Soc. Sci. Briixelles 99, 49 (1985). 

[91] M.J. Buckingham, in Phase Transitions and Critical Phemomena Vol. 2, edited by C. Domb, M.S. 
Green (Academic, London, 1972). 

[92] V.B. Magalinsky, Ya.P. Terletsky, Application of the microcanonical distribution to the statistical 
theory of multiple production of particles, Zh. Eksp. Teor. Fiz. 32, 584 (1957) [Sov. Phys. JETP 5, 483 
(1957)]. 

[93] J.I. Kapusta, Ch. Gale, Finite-Temperature Field Theory (Cambridge University Press, Cambridge, 
2006), Appendix. 

[94] J. A. McLennan, Jr., Statistical mechanics of the steady state, Phys. Rev. 115, 1405 (1959). 
[95] B. Robertson, Equations of motion in nonequilibrium statistical mechanics, Phys. Rev. 144, 151 (1966). 
[96] R.A. Treumann, C.H. Jaroschek, Gibbsian theory of power-law distributions, Phys. Rev. Lett. 100, 
155005 (2008). 



22 



